eLIFE 

elifesciences.org 



RESEARCH ARTICLE 



9 © 



Systems analysis of the CO2 concentrating 
mechanism in cyanobacteria 

Niall M Mangan^*, Michael P Brenner* 

School of Engineering and Applied Sciences and Kavli Institute for Bionano Science 
and Technology, Harvard University, Cambridge, United States 



Abstract Cyanobacteria are photosynthetic bacteria with a unique CO2 concentrating mechanism 
(CCM), enhancing carbon fixation. Understanding the CCM requires a systems level perspective 
of how molecular components work together to enhance CO2 fixation. We present a mathematical 
model of the cyanobacterial CCM, giving the parameter regime (expression levels, catalytic rates, 
permeability of carboxysome shell) for efficient carbon fixation. Efficiency requires saturating 
the RuBisCO reaction, staying below saturation for carbonic anhydrase, and avoiding wasteful 
oxygenation reactions. We find selectivity at the carboxysome shell is not necessary; there is an 
optimal non-specific carboxysome shell permeability. We compare the efficacy of facilitated 
CO2 uptake, CO2 scavenging, and HCO3" transport with varying external pH. At the optimal 
carboxysome permeability, contributions from CO2 scavenging at the cell membrane are small. We 
examine the cumulative benefits of CCM spatial organization strategies: enzyme co-localization and 
compartmentalization. 
DOI: 10.7554/eLife.02043.001 
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Introduction 

Intracellular compartments are critical for directing and protecting biochemical reactions. One of the 
simplest and most striking known examples of compartmentalization are the carboxysomes (Cannon 
et a/., 2001; Yeates et al., 2008) of cyanobacteria and other autotrophic proteobacteria {Savage et a/., 
2010; Rosgaard eta/., 2012). These small, 100-200 nm compartments separate the principal reaction 
of the Calvin cycle, the RuBisCO catalyzed fixation of carbon dioxide (CO2) into 3-phosphoglycerate, 
from the rest of the cell {Cannon et al., 1991). CO2 and oxygen (O2) competitively bind as substrates 
of RuBisCO, and the reaction with O2 produces phosphoglycolate, a waste product which must be 
recycled by the cell {Jordan and Ogren, 1981; Tcherkez et al., 2006; Savir et al., 2010). To maximize 
carboxylation and minimize oxygenation, the carboxysome is believed to act as a diffusion barrier to 
CO2 {Reinhold et al., 1989; Dou et al., 2008). There is much interest in the design and function of 
such compartments and whether they can be used to enhance carbon fixation in other organisms such 
as plants or to improve reaction rates in other metabolic systems {Papapostolou and Homorka, 2009; 
Agapakis et al., 2012; Ducat and Silver, 2012; Frank et al., 2013). Increased efficiency of biochem- 
ical reactions will lead to better yield in bioengineered bacterial systems, creating new possibilities for 
production of high-value products such as biofuels. Enhancing carbon fixation in plants or other organ- 
isms could lead to increased carbon sequestration, or crop yield. 

The concentrating mechanism in cyanobacteria relies on the interaction of a number of well charac- 
terized components, as shown In Figure 1, transferring inorganic carbon from outside the cell Into 
cytosol and carboxysomes {Allen, 1984; Kaplan and Reinhold, 1999; Badger and Price, 2003; 
Price et al., 2007). Due to this mechanism, inorganic carbon concentration is elevated well above 
200-300 jjM, the CO2 concentration required for saturating the RuBisCO. Additionally a high CO2 
concentration increases the ratio of CO2 to O2 so that carboxylation dominates over oxygenation. 
Concentrations of 20-40 mM inorganic carbon, up to 4000-fold higher than external levels, have been 
observed {Sultemeyer et al., 1995; Price et al., 1998; Kaplan and Reinhold, 1999; Woodger et al.. 
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eLif e digest Cyanobacterla are microorganisms that live in water and, like plants, they capture 
energy from the sun to convert carbon dioxide into sugars and other useful compounds. This 
process — called photosynthesis — releases oxygen as a by-product. Cyanobacterla were crucial in 
making the atmosphere of the early Earth habitable for other organisms, and they created the vast 
carbon-rich deposits that now supply us with fossil fuels. Modern cyanobacterla continue to sustain 
life on Earth by providing oxygen and food for other organisms, and researchers are trying to 
bioengineer cyanobacterla to produce alternative, cleaner, fuels. 

Understanding how cyanobacterla can be as efficient as possible at harnessing sunlight to 'fix' 
carbon dioxide into carbon-rich molecules is an important step in this endeavor. Carbon dioxide 
can readily pass through cell membranes, so instead cyanobacterla accumulate molecules of 
bicarbonate inside their cells. This molecule is then converted back into carbon dioxide by an 
enzyme found in specials compartments within cells called carboxysomes. The enzyme that fixes the 
carbon is also found in the carboxysomes. However, several important details in this process are not 
fully understood. 

Here, Mangan and Brenner further extend a mathematical model of the mechanism that 
cyanobacterla use to concentrate carbon dioxide in order to explore the factors that optimize 
carbon fixation by these microorganisms. Carbon fixation is deemed efficient when there is more 
carbon dioxide in the carboxysome than the carbon-fixing enzyme can immediately use (which also 
avoids wasteful side-reactions that use oxygen instead of carbon dioxide). However, there should 
not be too much bicarbonate, otherwise the enzyme that converts it to carbon dioxide is overwhelmed 
and cannot take advantage of the extra bicarbonate. 

Mangan and Brenner's model based the rates that carbon dioxide and bicarbonate could move 
in and out of the cell, and the rates that the two enzymes work, on previously published experiments. 
The model varied the location of the enzymes (either free in the cell or inside a carboxysome), and 
the rate at which carbon dioxide and bicarbonate could diffuse in and out of the carboxysome (the 
carboxysome's permeability). Mangan and Brenner found that containing the enzymes within a 
carboxysome increased the concentration of carbon dioxide inside the cell by an order of magnitude. 
The model also revealed the optimal permeability for the carboxysome outer-shell that would 
maximize carbon fixation. 

In addition to being of interest to researchers working on biofuels, if the model can be adapted 
to work for plant photosynthesis, it may help efforts to boost crop production to feed the world's 
growing population. 
DOI: 10.7554/eLife.02043.002 



200S). The bilipid outer and cell membranes are highly permeable to small uncharged molecules such 
as CO2 {Gutknecht et al., 1977; Missner et al., 2008), so instead the cell primarily accumulates the 
charged and less membrane soluble bicarbonate (HCO3 ) {Volokita et al., 1984; Price and Badger, 
1989). Active transporters, both constitutive and inducible, bring HCO3" into the cell {Omata et al., 
1999; Price et al., 2004, 2008), and mechanisms exist to actively convert CO2 to HCO3" at the thyla- 
koid and cell membrane {Shibata et al., 2001; Maeda et al., 2002; Price et al., 2008). Once it pas- 
sively diffuses into the carboxysome, HCO3" is rapidly brought into equilibrium with CO2 by the enzyme 
carbonic anhydrase, resulting in the production of CO2 near RuBisCO. Carbonic anhydrase is known to 
be localized on the interior side of the carboxysome shell {Cannon etal., 1991; Cot et al., 2008; Long 
et al., 2008; Yeates et al., 2008). The carboxysome shell must be permeable enough to allow HCO3" 
and 3-phosphoglycerate to readily diffuse in and out. The function of this system and its ability to 
concentrate inorganic carbon depends on the interplay between these various molecular components. 
Without a model, flux measurements cannot determine the components relative roles in enhancing 
the CO2 concentration in the carboxysome. To date, it has not been possible to directly measure the 
partitioning of the internal carbon concentration in the cytosol and carboxysomes. We wish to charac- 
terize the distribution of internal carbon. Visualizations of the location of the carboxysomes with fluo- 
rescent microscopy in S. elongatus PCC7942 demonstrated that the carboxysomes are evenly spaced 
along the centerline of the cell, {Savage et al., 2010), raising the question of how spatial organization, 
beyond simple partitioning, changes the efficacy of the system. 
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Figure 1. Schematic of the CCM in cyanobacteria. Outer and cell mennbranes (in black), as well as, thylakoid 
membranes where the light reactions take place (in green) are treated together Carboxysomes are shown as four 
hexagons evenly spaced along the centerline of the cell. The model treats a spherically symmetric cell, with one 
carboxysome at the center Active HCO3" transport into the cell is indicated (in light blue), as well as active 
conversion from CO2 to HCOa" at the membranes, sometimes called 'facilitated uptake' or 'scavenging' (in orange). 
Both CO2 and HCO3" can leak in and out of the cell, with CO2 leaking out much more readily. Both species 
passively diffuse across the carboxysome shell. Carbonic anhydrase (red) and RuBisCO (blue) are contained in the 
carboxysomes and facilitate reactions as shown. 
DOI: 10.7554/eLife.02043.003 



The goal of this study is to further develop a mathematical model of the CCM {Reinhold et a/., 
1989, 1991; Fridlyand et a/., 1996) that uses recent experimental progress on the CCM to untangle 
the relative roles of the different molecular components, and predict the region of parameter space 
where efficient carbon fixation occurs. We are considering conditions where CO2 is limiting (15 |jM 
external inorganic carbon) and, for the moment, ignore other biological pressures. In this context, 
efficient carbon fixation requires two conditions: First, the CO2 concentration must be high enough 
that RuBisCO is saturated, and the competitive reaction with O2 is negligible. We emphasize that for 
the oxygenation reaction to be negligible the CO2 concentration should be higher than needed to 
merely saturate RuBisCO. Secondly, the carbonic anhydrase within the carboxysome must be unsatu- 
rated, so that extra energy isn't wasted transporting unused HCO3" into the cell. 

Examination of the system performance with varying expression levels of HCO3" transporters, car- 
boxysome permeability, and conversion from CO2 to HCO3 , reveals a parameter window where these 
conditions are simultaneously satisfied. We comment on the relation of this window to measured car- 
bon pools, carbon fixation rates, and HCO3" transporters. We find that the HCO3" concentration in the 
cytosol is constant across the cell, set by the HCO3" transport and leakage rates, and depends very 
little on the carboxysome permeability. The carboxysome permeability does, however, set how the 
CO2 is partitioned between the carboxysome and cytosol. At optimal carboxysome permeability, 
HCO3" diffusion into the carboxysome is fast enough to supply inorganic carbon for fixation, but the 
rate of CO2 leakage out of the carboxysome is low. We explore the fluxes from CO2 facilitated uptake 
and scavenging with varying ratios of external CO2 and HCO3 . Finally we discuss the proportion the 
carbon concentration comes from different methods of spatial organization such as co-localization, 
encapsulation, and spatial location of carboxysomes. Concentration of carbonic anhydrase increases 
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the maximum rate of reaction for carbonic anhydrase per volume, causing carbonic anhydrase to sat- 
urate at a higher level of HCO3 , and achieve an order of magnitude higher local CO2 concentrations. 
Encapsulation of the reactions in an optimally permeable carboxysome shell achieves another order of 
magnitude of CO2 concentration. 

Reaction diffusion model 

We present our mathematical model, which captures all aspects of the CCM as described above. This 
model is an expansion of previously developed models {Reinhold et al., 1989, 1991; Fridlyand et a/., 
1996). Our three dimensional model of the CCM solves for both the HCO3" and CO2 concentration 
throughout the a spherical cell. We solve this model numerically and analytically at steady state for 
three different spatial organizations of carbonic anhydrase and RuBisCO in the cell (Figure 6): enzymes 
distributed evenly throughout the cell, enzymes localized to the center of the cell but not encapsu- 
lated (as they would be on a scaffold), enzymes encapsulated in a carboxysome. We compare the 
effects of these scenarios in the discussion section, and for now consider a spherical cell of radius 
Rb = 0.5 [jm with a single spherical carboxysome of radius = 50 nm containing RuBisCO and car- 
bonic anhydrase. Numerical computations are carried out with finite difference methods in MATLAB. 
The details of analytic solutions are given in the Supplementary file 1. 

We include the effects of 3D diffusion, active transport and leakage through the cell membrane, 
and reactions with carbonic anhydrase and RuBisCO. In the carboxysome (r < RJ, the equations gov- 
erning the HCO3" and CO2, here H and C respectively, are 

dtC=Dy^c+RcA-RRub (1) 



dtH= DV^H-RcA, (2) 

where here D is the diffusion constant, and Rca is the carbonic anhydrase reaction, and R^ub is the 
RuBisCO reaction. The carbonic anhydrase reactions follows reversible Michaelis-Menten kinetics 
{Kaplan and Reinhold, 1999; Price et al., 2007), 

f^CA (rl' W 

Kba^ca + KcaH + KbgC (3) 

where Vca and Vta are hydration and dehydration rates, proportional to the local carbonic anhydrase 
concentration. /<„ and are the concentration at which hydration and dehydration are half max- 
imum. The RuBisCO reaction follows Michaelis-Menten kinetics with competitive binding with O2, 
R/jufa = VmaxCI (C + Km), where = (1 -1- O/K.). Here V™,, is the maximum rate of carbon fixation and 
Km is the apparent half maximum concentration value, which has been modified to include competitive 
binding with O2, O. K is the dissociation constant of O2 with the RuBisCO and is the half maximum 
concentration with no O2 present. RuBisCO also requires ribulose-1 ,5-bisphosphate, the substrate 
which CO2 reacts with to produce 3-phosphoglycolate. Under CO2 limiting conditions it has been 
shown that there is sufficient ribulose-1 ,5-bisphosphate to saturate all RuBisCO active sites, and 
the reaction rates are independent of ribulose-1 ,5-bisphosphate concentrations {Mayo et al., 1989; 
Whitehead et al., 2014). 

In the cytosol there is no carbonic anhydrase or RuBisCO activity, so Rqa = 0 and RR„b = 0, and there 
is only diffusion of CO2 and HCO3 . We do not include the natural, but slow, interconversion of CO2 
and HCO3" in the cytosol. This assumption is a good one given that the HCO3" concentration is known 
to be held out of equilibrium in the cell {Volokita et al., 1984; Price and Badger, 1989). In agreement 
with this experimental observation, we find that all the other processes effecting the concentration of 
HCO3" in the cytosol happen much faster than the natural interconversion. 

Boundary conditions proscribe the inorganic carbon fluxes into the cell and the diffusion across the 
carboxysome boundary. We treat the inorganic carbon fluxes at cell and thylakoid membranes 
together. At this cell boundary, there is passive leakage of both CO2 and HCO3 : the velocity of CO2 
across the cell membrane, is about 1000-fold higher than that of HCO3 , k'^, due to the lower per- 
meability of the membrane to charged molecules. To account for active import of HCO3 , we combine 
the total HCO3" flux, from all HCO3" transporters. These transporters include BCT1 (encoded 
by cpm), which is thought to be powered by ATP; and BicA and SbtA which are thought to be 
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symporters between HCO3" and Na*, driven by the highly controlled electrochemical gradient for 
Na+ {Omata et a/., 1999; Price et al., 2004, 2008). Additionally, there are two complexes NDH-I3 and 
NDH-I4 responsible for converting CO2 to HCO3 . This conversion is thought to either decrease CO2, 
creating a gradient across the membranes and 'facilitating uptake' of CO2, or 'scavenge' CO2 which 
has escaped from the carboxysome. These are localized to the thylakoid and possibly the plasma 
membrane. They have been linked to the photosynthetic linear and cyclic electron transport chain 
{Shibata et al., 2001; Maeda et al., 2002; Price et al., 2008). It has been proposed that electron 
transport drives the formation of local alkaline pockets where CO2 more rapidly converts to HCO3 . We 
simply describe the conversion with a maximal reaction rate a, and concentration of half maximal 
activity of Ka. Combining these effects, the boundary condition setting diffusive flux of HCO3" and CO2 
at the cell membrane is 

DdrC = - /^y-' -H /C^ (Co., - Q,,o3o/) 

^cytosol M) 



DdrH - jcHout + — ^ + km (Hout ~ Hcytosol) 

where the subscript cytosol and out indicate we are taking the concentration immediately inside and 
outside the cell boundary respectively. The diffusion constant times partial derivatives with respect to 
the radial coordinate, r, define the diffusive flux at the membrane. 

The carboxysome shell is composed of proteins with a radius = 50 nm. As of yet, there have been 
no direct measurements of the carboxysome permeability to small molecules. Using the carboxysome 
geometry, we can calculate an upper bound for the diffusive velocity across the carboxysome shell, 
which is directly related to the carboxysome permeability. Crystal structures {Yeates et al., 2007, 
2008; Cheng et al., 2008) show the surface has approximately Npores = 4800 small pores with 
radius rpore~0.35 nm, and thickness / = 1.8 nm. If is the characteristic velocity that small molecules 
pass through the shell, these numbers imply the upper bound for diffusive transport 

k < — (N ) = 0.02-^. This calculation is done by taking the probability that a molecule will 



encounter a pore on the carboxysome shel 



N 3 X poresurfacearea 



and multipling it by the speed 



carboxysome surface area 

a small molecule will diffuse through the length of the pore (D//). Since it does not take into account 
any charge effects, which would add an additional energy barrier, it is an upper bound. Although there 
has been much speculation that the positively charged pores might enhance diffusion of negatively 
charged HCO3 [Cheng et al., 2008; Dou et al., 2008; Yeates et al., 2008), here we explore the sim- 
plest assumption, that both HCO3" and CO2 have the same permeability. Namely, the boundary condi- 
tions at the carboxysome shell are 



Dd^C kc {Ccytosol ^carboxysome) 



(6) 



DdrH ^c(^cytoso/ H i-gi-boxysome) • 



(7) 



We will vary (henceforth called carboxysome permeability, although it is a velocity) within our 
model and see that there is a range of Icc where the CCM is effective even with k^ identical for CO2 and 
HCO3-. 



Results 

Analysis of model: finding functional parameter space 

Now that we have defined our model, we wish to find the range of parameters where efficient carbon 
fixation occurs. In what follows, we fix the enzymatic rates, cell membrane permeability, and diffusion 
constant as reported in the literature {Gutknecht et al., 1977; Jordan and Ogren, 1981; Heinhorst 
et al., 2006; Missner et al., 2008; Tables 1 and 2). Note that full analytic solutions are available in 
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Table 1. 


Parameter values chosen for main set of simulations. 


unless otherwise 


indicated 


Parameter Definition 


Value 


Reference 


^out 


concentration of bicarbonate outside the cell 


14 pM* 


[Price et al., 2008) 


Cout 


concentration of carbon dioxide outside of cell 


0.14 pM* 


[Price et al., 2008) 


D 


diffusion constant of small molecules, CO2 and HCO3 


10-^ 
s 


[Fridlyand et al.. 1996) 




permeability of cell membrane to CO2 


cm 

0 3 


[Gutknecht et al., 1977, 
Missner et al., 2008) 


k" 


permeability of cell membrane to CO2 


3x10^ — 
s 


(Gutknecht et at 1 977- 
Missner et al., 2008) 


Rc 


radius of carboxysome 


5x1 0"' cm 


[Schmid et al., 2006, 
Cheng et al., 2008) 


D 

"b 


IdLllUb Ul LJdCLcrlld 


5x10"^ cm 


[Savage et al., 2010) 


ic 


HCO3" transport rate resulting in 30 mM cytosolic 
HCO3- pool 


0.7—* 
s 


calcu lated 


K 


optimal carboxysome permeability 


6x10-^* 
s 


calculated 




cell volume 


5.2 X 10-'° 


calculated 




cell surface area 


3 X 10-8 crr,2 


calculated 


*these parameters are varied in the text, but these values are used unless noted otherwise. 
DOI: 10.7554/eLife.02043.004 



Supplementary file 1 sections S3 and S4, so the effect of varying these parameters can be analyzed. 
We consider the efficacy of the CCM as a function of j„ the flux of HCO3" into the cell, k„ the carboxy- 
some permeability, and the parameters (a, KJ governing the CO2 facilitated uptake mechanism. Both 
a and can be regulated by the organism and vary depending on environmental conditions, whereas 
the carboxysome permeability, k^, is the parameter with the largest uncertainty and debate {Cannon 
et al., 2001, Cheng et al., 2008; Yeates et al., 2008). 

For any given pair of and j„ we ask whether the CO2 concentrating mechanism is effective, using 
the criteria of saturating RuBisCO, reducing oxidation reactions, and not increasing the HCO3" concen- 
tration beyond carbonic anhydrase saturation. Our central result is presented in Figure 2, which shows 
the range of k^ and where these conditions are met, assuming that there is no facilitated uptake, 
a = O.The blue shaded region shows where RuBisCO is unsaturated, and the red shaded region 
shows where carbonic anhydrase is saturated. There is a crescent shaped region between these 
regions, where the CCM is effective according to our criteria. In the white region oxygenation 



Table 2. Table comparing enzymatic rates {Sultemeyer et al., 1995; Woodger et al., 2005; 
Heinhorst et al., 2006) 



Enzyme reaction 


active 
sites 


^cut 


1 

s 


V™, in cell 


s 


V„a„ in carboxysome 


s 


Kv2 [hM] 


carbonic anhydrase 
hydration 


80 


8 X 10* 




1.1x10" 




2 x 10" 




3.2 X 103 


carbonic anhydrase 
dehydration 


80 


4.6 X 10* 


9.5x10" 




1.7 X 10" 




9.3 X 10^ 


RuBisCO 
carboxylation 


2160 


26 






103 




1.7 X 10' 




270 


in cell and carboxysome refer to the volumetric reaction rate assuming the enzymes are distributed through 
cellout the entire cell or only carboxysome. Vt, [V,^,^ for carbonic anhydrase dehydration) is estimated by assuming 
Keq = 5 and using parameters found in Heinhorst et al. (2006). V^, is V„,,„ for carbonic anhydrase hydration. 
Similarly, Kts, and /<„ are /C,/2 for dehydration and hydration respectively. 



DOI: 10.7554/eLife.02043.005 
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CO 

E 
o 



10' 



CC 10 
CD 

E 
i_ 

(D 
Q. 

(U 

E 
o 

CO 

>^ 

X 

o 

CO 

O 



10' 



1 


I Carbonic anhydrase 
I saturated 


Co 


I HCO = 30 mM in cell 

3 


RuBisCON^ 
unsaturated 


< 1 % error ^^V^ 
0^ reactions 



10"' 



10' 



10" 



10^ 



10' 



10^ 



Active HCO^" transport (cm/s) 



Figure 2. Phase space for HCO3" transport, and carboxysome permeability Plotted are the parameter values 
at which the CO2 concentration reaches some critical value. The left most line (dark blue) indicates for what values 
of jc and the CO2 concentration in the carboxysome would half-saturate RuBisCO (KJ. The middle line (light 
blue) indicates the parameter values which would result in a CO2 concentration where 99% of all RuBisCO reactions 
are carboxylation reactions and only 1% are oxygenation reactions when O2 concentration is 260 [iM. The right 
most (red) line indicates the parameter values which result in carbonic anyhdrase saturating. Here a = 0, so there is 
no CO2 scavenging or facilitated uptake. The dotted line (grey) shows the k^ and values, where the HCO3" con- 
centration in the cytosol is 30 mM. The HCO3" concentration in the cytosol does not vary appreciably with in this 

cm 

parameter regime, and reaches 30 mM at = 0.7 — . All other parameters, such as reaction rates are held fixed 
and the value can be found in the Table 1 and 2. 
DOI: 10.7554/eLife.02043.006 

The following figure supplements are available for figure 2: 

Figure supplement 1. Phase space for HCOs'transport and carboxysome permeability. 
DOI: 10.7554/eLife.02043.007 

Figure supplement 2. Effect of C02 scavenging or facilitated uptake on phase space for HCOs'transport, and 
carboxysome permeability, K^. 
DOI: 10.7554/eLife.02043.008 



reactions happen at a rate of greater than 1%. In the green shaded region oxygenation reactions 
occur at a rate of less than 1%. Within the white and green regions the CO2 concentration in the car- 
boxysome varies greatly. 

The lines dividing the regions in Figure 2 are lines of constant carboxysomal CO2 concentration in 
jc and Icc parameter space. The dark blue line is where CO2 = K^, the CO2 concentration for half max- 
imum RuBisCO reactions. The light blue line indicates parameter values resulting in the CO2 concen- 
tration (C99%) where rate of oxygenation reactions is 1% for O2 concentration of 260 Varying 
carboxysome permeability, values, require more or less HCO3 transport, to achieve the same 
carboxysomal CO2 concentration. 

We can calculate an amplification factor for the Cm line of constant carboxysomal CO2 concen- 

tration as Ac = "'"^"''y^"'"^ =133. Any combination of jcand Icc which produce C= Gn%, make 133 times 

more CO2 available in the carboxysome than there is total inorganic carbon outside the cell. Generally, 
increasing HCO3" transport, below the carbonic anhydrase saturation point results in higher CO2 con- 
centration in the carboxysome. 
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Varying HCO3" transport saturates 
enzymes 

The basic physics of the phase diagram Figure 2 
follows from examining how CO2 and HCO3" in the 
carboxysome change as is varied. Figure 3 shows 
the response to varying j„ with kc = 6 x 10" — 
(the optimal value in Figure 4). ^ 
When is low, the ratio of CO2 and HCO3" is 
constant, set by the chemical equilibrium at a given 
pH. In this case the rate of the carbonic anhydrase 
reaction is much faster than diffusion within the 
carboxysome, so that VbaKcaH= VcaKtaC. Unlike 
the uncatalyzed interconversion of CO2 and 
HCO3" in the cytosol, carbonic anhydrase brings 
the concentrations in the carboxysome to equilib- 
rium very quickly. The chemical equilibrium is 
Keq = H/C = (KfaaVca) (KcaVfaa) = 5, for pH around 7 
{DeVoe and Kistiakowsky, 1961), so that 
HCO3 >C02 in the carboxysome. Increased pH 
would increase and the proportion of 

HCO3 , while decreased pH would decrease K^^ 
and the proportion of HCO3 . Such variations 
do not substantially effect the subsequent dis- 
cussion and mechanisms, although they will 
change the absolute values of CO2 concentration 
in the carboxysome. 

The dashed line in Figure 3 shows the CO2 
level above which RuBisCO reaction is saturated: 
this gives the RuBisCO saturated (blue) boundary 
in Figure 2. We have similarly marked the con- 
centration €<}<)% where there is a 1% oxygen reac- 
tion error rate with a dash-doted line. 

At higher levels, the CO2 concentration no 
longer increases with increasing j„ because the 
carbonic anhydrase is saturated. The saturated 
regime occurs in Figure 3 when H„rbo)iysome> ^ba, so that increasing H^arboxysome (controlled directly by jj 
no longer increases the rate of production of Ccarboxysome- This transition from unsaturated to saturated 
carbonic anhydrase defines the line for the carbonic anhydrase saturated region in Figure 2. 

Carboxysome permeability has optimal value 

For each line of constant concentration in Figure 2 there is an optimal permeability value, where the 
least HCO3" transport is required to achieve the same CO2 concentration. The optimal permeability 
value shifts downward with increasing CO2 concentration (compare light and dark blue curves). For 

€99% the optimal permeability is kc = 6 x 10~'^— , below the calculated upper bound: kc < 0.02 — 

s s 

obtained above from the carboxysome structure. To further understand the effect of permeability, we 
examine the CO2 concentration in the carboxysome for varying carboxysome permeabilities and a 
fixed HCO3" transport rate in Figure 4. Figure 4A, shows that there is a broad range of where the 
CCM has maximal efficacy. Figure 4 shows the distribution of inorganic carbon throughout the cell 
when the permeability is low (B), optimal (C), and high (D). At high permeability, the CO2 produced in 
the carboxysome rapidly leaks out of the carboxysome, and the CO2 concentration in the cytosol, 
shown in Figure 4D, is relatively high. When the carboxysome permeability decreases to near the 
optimal value. Figure 4D, the carboxysome traps CO2, and the cytosolic levels are lower, decreasing 
leakage out of the cell. This transition occurs when diffusion across the cell (and carboxysome) takes a 
shorter time than diffusion through the carboxysome shell; or the CO2 in the carboxysome is effec- 
tively partitioned from the CO2 in the cell. 
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Figure 3. Numerical solution (diamonds and circles) 
and analytic solutions (carbonic anhydrase unsaturated, 
solid lines, and saturated, dashed lines) correspond 
well. HCO3" transport is varied, and all other system 
parameters are held constant. The CO2 concentration 
above which RuBisCO is saturated is K„ (grey dashed 
line). The CO2 concentration where the oxygen 
reaction error rate will be 1% is Cm. (grey dash-dotted 
line). The transition between carbonic anyhdrase 
being unsatruated and saturated happens where the 
two analytic solutions meet (where the dashed and 
solid red lines meet). Below a critical value of transport, 

]^ le"^-^ the level of transport is lower than the 
s 

HCO3" leaking through the cell membrane. 
DOI: 10.7554/eLife.02043.009 

The following figure supplements are available for 
figure 3: 

Figure supplement 1. No effect of localizing carbonic 
anhydrase to the shell of the carboxysome. 
DOI: 10.7554/eLife.02043.010 
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Figure 4. Concentration of CO2 in the carboxysome with varying carboxysome permeability. (A) Numerical solution 
(diamonds and circles) and analytic solutions (carbonic anhydrase unsaturated, solid lines, and saturated, dashed 
lines) correspond well. On all plots CO2 (red circle) < HCO3" (blue diamond). Concentration in the cell along the 

radius, r, with carboxysome permeability = ^e~^^'^ (B), = ^e~^^^ (C), k^ =1"^ (D). Grey dotted lines in (B), (C), 

s s s 

(D) indicate location of the carboxysome shell boundary. The transition from low CO2 at high permeability (D) to 

maximum CO2 concentration at optimal permeability (C) occurs at k* = — = 2-^^. At low carboxysome 

R, s 

cm 

permeability (B) HCO3" diffusion into the carboxysome is slower than consumption. For all subplots a = 0 and 

s 

/ = 0.7 . Qualitative results remain the same with varying increasing a will increase the gradient of CO2 across 

s 

the cell as CO2 is converted to HCO3" at the cell membrane. 
DOI: 10.7554/eLife.02043.011 



If the carboxysome permeability is below optimal, diffusion of HCO3" Into the carboxysome cannot 
keep up with consumption from RuBisCO. The existence of an optima requires RuBisCO consumption 
to be low enough that there is a where the cytosol and carboxysome are partitioned, but HCO3" 
diffusion in can keep up. When such an optima exists, the carboxysome permeability can improve 
the CO2 concentration in the carboxysome without any special selectivity between HCO3" and CO2. 
The location and concentrating power of the optimal regime, is dependent on the size of the cell 
and the membrane permeabilities to CO2 and HCO3 . 

Discussion 

Are the fluxes and concentrations reasonable? 

While we have solved our model to describe a vast parameters space it is instructive to compare the 
fluxes and concentrations we find within our optimal parameter space (the green region in Figure 2) 
to actual numbers. At low external inorganic carbon conditions, internal inorganic carbon pools due to 
CCM activity are regularly measured as high as C, = 30mM. The inorganic carbon is predominantly in 
the form of HCO3 , and measurements do not distinguish between the cytosol and carboxysome 
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{Sultemeyer et a/., 1995; Price et a\., 1998, 2008; Kaplan and Reinhold, 1999; Woodger et a/., 

cm 

200S). In our model, we find that the cytosolic HCO3" concentration is 30 mM when jc - 0.7 — , over 

s 

a wide range of the carboxysome permeability (indicated as the dashed grey line in Figure 2). 

From Figure 4 we can also see that the cytosolic HCO3" concentration is the dominate form of inor- 
ganic carbon in the cell at jc = 0.7 — . We examine the fate of the HCO3" transported into the cell in 

s 

terms of the HCO3" leaking out, CO2 leaking out, CO2 fixation or carboxylation, and O2 fixation or 
oxygenation {Table 3). 

For cells grown under low inorganic carbon conditions net HCO3" fluxes (transport-leakage) are 

measured 10^ prr\o\ ^ qq^ ^^^^ being slightly lower but the same order of magnitude 

(mgChI s) 

{Badger et al., 1994; Whitehead et al., 2014). High external inorganic carbon conditions produce 

slightly higher net HCO3" rates {Tchernov et al., 1997). Assuming chlorophyll per cell volume of 

around 10~^^ "^^^^^ for cells of our size we can convert this into a flux of 10~^ pf^ol ^^^^^^^ 
cell (cells) 

2002, 2004; Whitehead et al., 2014). The net HCO3- flux for our model cell is 10"^ pmol 

(cells) 

about an order of magnitude too high. If we choose a HCO3" transport value one order of magni- 
tude smaller, we will get net fluxes of the same order of magnitude as the measurements at the 
cost of slightly lower carboxylation rates and higher oxygenation rates {Table 4). This would also 
mean a lower internal HCO3" pool. Alternatively, the same internal HCO3" could be reached at a 
lower flux rate, if the external HCO3" is higher. We have chosen a dramatically low external inorganic 
carbon concentration, where the CCM is known to be up-regulated (Price et al., 2008). The general 
results we present hold until the internal and external concentration are of the same order of magni- 
tude, at which point the CCM is no longer necessary. Since the majority of the HCO3" transport is 
balanced by HCO3" leakage, we can find the transport rate needed to sustain a particular internal 
concentration by the simple formula: jc = /Cm(Hout~ Hcytosoi(Rb))/Hout. 

While we can compare the net fluxes, we have not found direct experimental evidence for the 
absolute HCO3" uptake rate. To determine whether this HCO3" transport rate is reasonable we 
perform a back of the envelope calculation. Our simulated cell has a flux of 2 x 1 0^molecules/s. 

1 r r ^^-^ molecules , ,,, r 1 • 

Assuming the rate of transport per transporter of 10 and our cells surface area this 

s 



Table 3. Fate of carbon brought into the cell for j^ = 0.7 cm/s and Ic^ = 6x10 ^cm/s 

pmol 



formula 



(cell s) 



% of flux 



HCO3 transport 



3.3 X 10- 



HCO3 leakage 



3.2 X 10-' 



96.6 



CO2 leakage 
carboxylation 

oxygenation 



C + K 



1 + 



V° C 



1 + - 



10-' 



9 X 10-' 



3 X 10-' 



3.1 

0.3 

0.001 



DOI: 10.7554/eLife.02043.012 
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Table 4. Fate of carbon brought into the cell for = 0.07 cm/s and =6x10"^ cm/s 

pmol 





formu 


a 




(cell s) 


% of flux 


HCO3 transport 








3.3 X 10-^ 




HCO3" leakage 








3.2 X 10-^ 


95.3 


CO2 leakage 




~Ccytoso/('^b)) 


9 X 10-' 


2.8 


carboxylation 








6 X 10-' 


1.7 














oxygenation 


v° 






2 X 10-8 


0.07 






[ Kc> 








DO!: 10.7554/eLife.02043.013 



requires about 10^ transporters ^^^.^ actually not that far off from the number of ATP synthase com- 
um I 

coitidIgxgs 

plexes on the thylakoid membrane in spinach, 700 ^ {Miller and Staehelin, 1979), 

although it is still quite high. 

According to our calculation only around 1 % of the inorganic carbon transported into the cell is 
fixed into 3-phosophoglycerate. The conclusion that about 99% of inorganic carbon transported into 
the cell is lost through leakage challenges the assumption that the 3 ATP and 2 NADPH used during 
the Calvin-Benson-Bassham cycle is the dominant energy expenditure. If it holds true, then cyanobac- 
teria invest much more energy in inorganic carbon uptake than was previously understood. Even in this 
highly CO2 concentrating regime, 5x10'' 2-phosophoglycolate are produced per second. Cyanobacteria 
have been shown to have multiple pathways for recycling 2-phosophoglycolate {Hackenberg et al., 
2009). Our system fixes CO2 at a rate of 0.14 pg/hr. Given the volume of our cell, and the fact that 
between 1 1 5-300 fg/jjm^ of carbon are needed to produce a new cyanobacterial cell {Mahlmann 
et al., 2008) we need between 0.1 and 0.3 picograms of carbon per cell. At the higher flux rate 
{Table 3) this means that a cell could replicate every 1-2 hr, so faster than cyanobacteria replicate. 
The lower flux rate {Table 4) would produce fix enough CO2 for the cell to replicate every 8 to 21 hr, 
which is similar to the division times of cyanobacteria. 

Concentration profiles of CO2 and HCOs" across the cell 

cm 

At jc = 0.7 — , varying the carboxysome permeability changes how the available inorganic carbon 
s 

is partitioned between the carboxysome and cytosol, thereby setting the carboxysomal CO2 concen- 
tration as shown in Figure 4. Strikingly, the HCO3" concentration is constant across the cytosol. This is 
because the cell membranes have low permeability to HCO3 ; thus, the rate of escape is slow and 
HCO3" equilibrates across the cell. A consequence of this flat HCO3" profile is that the carboxysome 
experiences the same HCO3" concentration, independent of its position in the cell. This means the 
incoming inorganic carbon source for the carboxysome system is invariant with the position of the 
carboxysome in the cell. 

In contrast, there is a gradient in CO2 concentration across the cell when the carboxysome perme- 
ability is at or above the optimum {Figure 4C,D). The cell membrane is more permeable to CO2. The 
gradient means that the CO2 leakage out of the cell affects the CO2 leakage out of the carboxysome. 
Moving the carboxysome close to the cell membrane increases the leakage rate of CO2 out of the 
carboxysome. Notably, in S. e/ongatus the carboxysomes are located along the center line of the cell, 
away from the cell membranes {Savage et al., 2010). The spatial profiles of HCO3" and CO2 give no 
hint as to why the carboxysomes are spaced apart from one another. Since the gradient in HCO3" is 
flat, there is no competition between the carboxysomes for HCO3" (the main incoming source of 
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inorganic carbon). In fact, since the local concentration of CO2 is higher near a carboxysome, nearby 
carboxysomes would "feed" each other CO2. As has been shown, such clumping would reduce the 
probability of distributing carboxysomes equitably to daughter cells, possibly counteracting any ben- 
efit {Savage et a/., 2010). 

The concentration in the carboxysome is basically constant, because the carboxysome is so small 
that diffusion across it takes very little time. A consequence of this is that the organization of the reac- 
tions in the carboxysome does not effect the CO2 concentration in the carboxysome {Figure 3 — figure 
supplement 1). Therefore, the localization of the carbonic anhydrase to the inner carboxysome shell 
seems to have no effect on the CCM. It has been suggested that diffusion in the carboxysome should 
be slower, since the carboxysome is packed with RuBisCO. One proposed consequence of slower 
diffusion in the carboxysome is that it could trap CO2, making a low carboxysome permeability unnec- 
essary. We have tested this hypothesis {Figure 2 — figure supplement 1), and find that assuming the 

diffusion constant one would expect for small molecules in a 60% sucrose solution {D^ = 10 '^^), 

s 

does reduce the optimal carboxysome permeability. However, for any carboxysome permeability a 
higher HCO3" transport rate is needed to achieve the same carboxsomal CO2 concentration. So, if the 
diffusion is slower in the carboxysome it does not aid the CCM. Even at this slower diffusion, the CO2 
concentration across the carboxysome is flat. 

Benefit of CO2 to HCOs" conversion: facilitated uptake or scavenging of 
CO2 

We investigate the effect of CO2 to HCO3" conversion at the thylakoid and cell membranes (combined 
in the model). Increasing conversion, a >0, can facilitate uptake of CO2 from outside the cell and 
scavenge CO2 escaped from the carboxysome. Facilitated uptake results in saturating both carbonic 
anhydrase and RuBisCO at a lower level of HCO3" transport. Scavenging broadens the range of car- 
boxysome permeability which will effectively separate the inorganic carbon pools in the carboxysome 
and outside. Scavenging decreases the concentration of CO2 in the cytosol, so a more permeable 
carboxysome can still result in a low leakage rate of inorganic carbon out of the cell (more of the 
inorganic carbon in the cytosol is in the form of HCO3" which leaks out less readily). However, neither 
of these effects is particularly strong in our current range of reaction rates, and cell membrane perme- 
ability {Figure 2 — figure supplement 2). 

The relative effects of these two mechanisms depends on the external CO2 and HCO3" concentra- 
tions. In saltwater environments the pH is near 8 and HCO3" is the predominant inorganic carbon 
source. While external pH is not explicitly treated in our model, we can account for changes to pH 
through the external inorganic carbon concentration. To be consistent with oceanic environment, thus 
far we have shown results for low external inorganic carbon concentrations of [CO2] = 0.1 \iM and 
[HCO3 ] = 14.9 |jM. The effect of facilitated uptake, under these assumptions, is very small. In freshwa- 
ter or under conditions of ocean acidification, where the pH could fall to 6 or lower, there can be a 
much larger proportion of CO2 (>50%). Figure 5 shows the absolute contribution of HCO3" transport, 
facilitated CO2 uptake, and CO2 scavenging for varying proportions of external CO2. Even though we 
assume the same velocity of facilitated uptake and HCO3" transport (j = — = 1^^), facilitated 

uptake contributes less because it is limited by CO2 diffusion across the membrane. At the same rates 
of transport the facilitated uptake mechanism only contributes more than active HCO3" if the CO2 
concentration is greater than 80% of external inorganic carbon. This is consistent with observations 
that oceanic cyanobacteria such as Prochlorococcus only seem to possess gene homologs for HCOs^ 
transport systems, while other freshwater and estuarine cyanobacteria have gene homologs for both 
constitutive (NDH-I4) and inducible (NDH-I3) CO2 uptake systems as well as inducible HCO3" trans- 
port systems (BicA, SbtA, and BCT1) {Price, 2011). 

Scavenging only contributes significantly to total incoming HCO3" when the carboxysome permea- 
bility is higher than optimal. Figure 5, and does not contribute significantly below our caculated upper 

bound of kc < 0.02—. In these ranges for carboxysome permeability, there is very little CO2 leaking 
out of the carboxysome into the cytosol, so there is very little CO2 to scavenge. Figure 5. The effect 
of scavenging is dependent on the cell membrane permeability to CO2 and HCO3 . 

Given that scavenging has no obvious affect on HCO3" concentrations, it is reasonable to wonder 
why this mechanism exists at all. One might assume that scavenging prevents leakage, but if the energy 
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to bring a 'new' HCO3" molecule from outside the 
cell is the same as the energy required to save an 
'old' CO2 molecule from leaking out, there is no 
obvious advantage of preventing the leakage. It is 
possible that since the scavenging mechanism is 
associated with the electron transport chain of the 
light reactions of photosynthesis scavenging can 
be ramped up more easily when there is excess 
light energy. If this were the case, a comparison of 

jc = 1 — and — = 1 — is deceiving and — could 

s Ka s Ka 
be much larger. It has been suggested that the 
cell uses this mechanism as a way to dissipate 
excess light energy {Tchernov et a/., 1997, 2003). 

Cellular organization 

The most striking aspect of the CCM is the way 
that spatial organization is used to increase the 
efficacy of the reactions. Figure 6 compares the 
effect of different enzymatic reaction organiza- 
tions. Concentrating carbonic anhydrase and 
RuBisCO to a small region in the center of the 
cell, on a scaffold for example, leads to an order 
of magnitude increase in the concentration of 
CO2. Localizing the carbonic anhydrase to a small 
volume concentrates it, increasing the max- 
imum reaction rate per volume, \/„ and Vt,a. A 
larger \4a increases the HCO3" concentration at 
which carbonic anhydrase is saturated allowing 
the mechanism to take advantage of a larger 
HCO3" flux, jc- A small increase can be gained 
from encapsulating the enzymes in a permeable 
carboxysome shell and another order of magnitude 
is gained at the optimal permeability. At optimal carboxysome permeability, the CO2 is effectively parti- 
tioned into the carboxysome and conversion can act only as facilitated uptake as shown in Figure 5. 

Another advantage of localizing the enzymes in a small region at the center of the cell is separating 
carbonic anhydrase from the a(C02to HCO3 ) conversion mechanism, preventing a futile cycle. The 
futile cycle is most detrimental when the enzymes are distributed throughout the cytosol, and increases 
the oxygenation error rate by 10% (data not shown). Concentrating the enzymes away from the cell and 
thylakoid membranes, where conversion happens, removes this effect. With a carboxysome or scaffold 
the oxygenation rate is almost exactly the same with and without the a conversion mechanism. This 
is consistent with the previously shown detrimental effect of having active carbonic anhydrase free 
within the cytosol {Price and Badger, 1989). It would be impossible to keep the cytosol completely free 
from carbonic anhydrase enzyme, so there must be a way of activating it within the carboxysome only. 
Carbonic anhydrase is inactivated under reducing conditions (Pena et a/., 2010). Recently it was 
shown that carboxysomes oxidize after assembly, providing a way to keep carbonic anhydrase inactive 
until fully enclosed in a carboxysome (Chen et a/., 2013). 

Effects of pH 

Cyanobacteria must regulate pH as almost all biochemical reactions are pH sensitive. We do not at- 
tempt to model this regulation or potential pH variation within the cell, however pH may be 
included implicitly in a couple ways. We have already explored the effect of varying external pH, 
and the effects of pH on carbonic anhydrase. Cytosolic pH would have little direct effect on the CO2 
and HCO3" levels since the non-enzymatic interconversion is very slow as previously discussed. The 
effect of internal pH could also be explored by varying the reaction rate of RuBisCO, which is pH sen- 
sitive. Varying the reaction rate of RuBisCO greatly could change the range where a non-specific 
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Figure 5. Size of the HCOa^ flux in one cell from varying 
sources, as the proportion of CO2 to HCO3" outside 
the cell changes changes. We show results for three 
carboxysome permeabilities, k^, and only the scavenging 

is effected. Total external inorganic carbon is 15 pM, 

.cm , a .cm , , 

/ =1 — and — = 1 . When the carboxysome 

s /<„ s 

permeability is larger than optimal, =1 , scavenging 

s 

can contribute more than facilitated uptake at low 
external CO2 concentrations. However, when the 
carboxysome permeability at or below our geometric 

, , k <0.02— . . ,. I II 

bound, s , scavenging IS negligibly small. 

Unless there is very little HCOa^ in the environment, 
HCOa" transport seems to be more efficient than CO2 
facilitated uptake. 
DOl: 10.7554/eLife.02043.014 
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carboxysome permeability can Increase the con- 
centration of CO2 in the carboxysome. It would 
be unexpected that the RuBisCO rate be much 
faster than we assume, as we have assumed a rate 
on the high end. A lower RuBisCO rate would 
increase the range of effective carboxysome per- 
meabilities. As previously mentioned it has been 
hypothesized that the CO2 facilitated uptake 
mechanism functions by creating local alkaline 
pockets. Diffusion of hydrogen ions across the cell 
would be very fast, so it could be very difficult to to 
maintain local alkalinity. Whether such pH gradients 
are possible, is certainly a subject of future interest. 

Conclusions 

We have described and analyzed a model for the 
CO2 concentrating mechanism in cyanobacteria. 
There exists a broad range of HCO3" transport 
and carboxysome permeability values which 
result in effective CO2 concentration in the car- 
boxysome. This effective concentration param- 
eter space is defined by CO2 levels high enough 
to saturate RuBisCO and produce a favorable 
ratio of carboxylation to oxygenation reactions, 
but not so high as to saturate carbonic anhy- 
drase (after which increasing HCO3" transport will 
not increase the CO2 concentration). An optimal, 
non-specific carboxysome permeability exists, 
where HCO3" diffusion into the carboxysome is 
not substantially inhibited, but CO2 leakage is 
minimal. HCO3" concentrations across the cell are flat and are predominately set by the transport rate in, 
and leakage out. We quantitatively compare the transport rates and concentrations we predict in our 
optimal parameter space, and find them to be in good agreement with experiment. We also comment 
on the effects of external pH on CO2 versus HCO3" uptake mechanisms. Finally we describe the cumu- 
lative benefits of co-localization, encapsulation, and optimal carboxysome permeability on the CCM. 

Further comparison of this model to experimental flux measurements, especially to determine the 
quantitative contributions of different transporters under different physiological conditions would be 
very interesting. Current solutions are for steady state at constant external concentration, but most 
gas exchange measurements, by necessity, measure the fluxes as the inorganic carbon is depleted in 
the media. The model could be modified to solve the time dependent problem with varying external 
inorganic carbon. As of yet the carboxysome permeability has not been measured directly, and it 
would be quite interesting to see how close it is to our 'optimal' prediction. 
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Figure 6. Concentration of CO2 achieved through various 
cellular organizations of enzymes, where we have selected 
the HCO3" transport level such that the HCOs^ concen- 
tration in the cytosol is 30 mM. O2 concentration is 260 ^iM. 
The oxygenation error rates, as a percent of total 
RuBisCO reactions are indicated on the concentration 
bars. The cellular organizations investigated are RuBisCO 
and carbonic anhydrase distributed throughout the 
entire cytosol, co-localizing RuBisCO and carbonic 
anhydrase on a scaffold at the center of the cell without 
a carboxysome shell, RuBisCO and carbonic anhydrase 
encapsulated in a carboxysome with high permeability 
at the center of the cell, and RuBisCO and carbonic 
anhydrase encapsulated in a carboxysome with optimal 
permeability at the center of the cell. 
DOI: 10.7554/eLife.02043.015 
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